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ABSTRACT 

Aims. To study the cosmological evolution of Active Galactic Nuclei (AGN) is one of the main goals of X-ray surveys. To accurately 
determine the intrinsic (before absorption) X-ray luminosity function, it is essential to constrain the evolutionary properties of AGN 
and therefore the history of the formation of supermassive black holes with cosmic time. 

Methods. In this paper we investigate the X-ray luminosity function of absorbed (log N H > 22) and unabsorbed AGN in three energy 
bands (Soft: 0.5-2 keV, Hard: 2-10 keV and Ultrahard: 4.5-7.5 keV). For the Hard and Ultrahard sources we have also studied the 
Nh function and the dependence of the fraction of absorbed AGN with luminosity and redshift. This investigation is carried out using 
the XMS survey along with other highly complete flux-limited deeper and shallower surveys in all three bands for a total of 1009, 
435 and 1 19 sources in the Soft, Hard and Ultrahard bands, respectively. We have modelled the instrinsic absorption of the Hard and 
Ultrahard sources (Nh function) and computed the X-ray luminosity function in all bands using two methods. The first makes use of 
a modified version of the classic l/V„ technique, while the second performs a Maximum Likelihood fit using an analytic model and 
all available sources without binning. 

Results. We find that the X-ray luminosity function (XLF) is best described by a Luminosity-dependent Density Evolution (LDDE) 
model. Our results show a good overall agreement with previous results in the Hard band, although with slightly weaker evolution. 
Our model in the Soft band present slight discrepancies with other works in this band, the shape of our present day XLF being 
significantly flatter. We find faster evolution in the AGN detected in the Ultrahard band than those in the Hard band. 
Conclusions. The results reported here show that the fraction of absorbed AGN in the Hard and Ultrahard bands is dependent on the 
X-ray luminosity. We find evidence of evolution of this fraction with redshift in the Hard band whereas in the Ultrahard band there is 
none, possibly due to the low statistics. Our best-fit XLF shows that the high-luminosity AGN, detected in all bands, exhibit a similar 
behaviours and are fully formed earlier than the less luminous AGN. The latter sources account for the vast majority of the accretion 
rate and mass density of the Universe, according to an anti-hierarchical black hole growth scenario. 
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1. Introduction 

One of the main goals of X-ray surveys is to study the cos- 
mological properties of Active Galactic Nuclei (AGN) as they 
are strongly linked to the accretion history of the Universe and 
the formation and growth of the supermassive black holes that 
are believed to reside in the centre of all galaxies, active or not 
(Kormendy & Richstone 19951 Magorrian et al. 1998 Richstone 
et al. I1998I I. The first studies in this field were constrained to 
the soft X-ray band (<2 keV) which could be biased against ab- 
sorbe d AGN (Maccacaro et al . [19911 Boyle et al. [TW3l Page 
et al. 119971 Miyaji et al. 12000ft . Therefore, hard X-ray surveys 
(>2 keV) are essential to describe the luminosity function of 
the AGN population, including obscured AGN which should 
be the main contributors to the cosmic X-ray background (see 
Setti & Woltjer [19891 and Fabian & Barcons [19921 for a review). 
Moreover, the large majority of energy density generated by ac- 
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cretion power seems to take place in obscured AGN (Fabian et 
al. U9981 l. as demonstrated by the integrated energy density of the 
cosmic X-ray background (see e.g. Comastri et al. 1 19951 Gilli et 
al. 120071 1. Ignoring the obscured AGN population could there- 
fore bias our understanding of the cosmic evolution of the struc- 
tures in the X-ray Universe. Furthermore, it has been confirmed 
that the fraction of absorbed AGN decreases with increasing X- 
ray luminosity in X-ray selected samples of AGN (Ueda et al. 
120031 Steffen et al. |20U3T > as well as in optically selected AGN 
(i.e. Simpson 2005). In addtion, some authors claim to have 
found a positive evolution of the fraction of absorbed AGN with 
redshift (La Franca et al. 120051 Ballantyne et al. 2006 Treister 
& Urry 2006 and, more recently, Hasinger 2008). The evolution 
of the fraction of the absorbed/obscured AGN is still a subject 
of debate and its study is therefore one of the primary targets of 
AGN surveys. 

Hard X-ray photons with energies between 2 and 10 keV can 
pass through large amounts of matter without being absorbed 
and hence are very useful to detect absorbed sources with in- 
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trinsic column densities up to Nh < 10 24 cirT 2 , although they 
are not energetic enough to escape from Compton-thick objects 
(Nh ^ 10 24 cm -2 ). Previous hard X-ray surveys in the 2-10 keV 
band (Cagnoni et al. 119981 Fiore et al. 119991 Giacconi et al. 
I200~n Baldi et al. [20021 Harrison et al. [20031 Alexander et al. 
2003) have provided least biased and highly complete samples 
of sources, many of them being obscured AGN. Notwithstanding 
this, an important fraction of the AGN detected in deep fields 
fail to provide good quality X-ray information, and their opti- 
cal counterparts are far too faint to make reliable spectroscopic 
identifications of the totality of the samples. 

Previous works have studied the X-ray luminosity function 
(XLF, hereafter) of AGN, mainly in the 0.5 -2 keV and 2-10 keV 
bands. For instance, Miyaji et al. (2000) and Hasinger et al. 
d2003T l studied the 0.5-2 keV XLF of Type- 1 AGN testing a vari- 
ety of models. Their results have ruled out both Pure Luminosity 
and Pure Density Evolution models in favour of a Luminosity 
Dependent Density Evolution model which best describes the 
observed XLF of these sources. In this paper we use a sample 
that contains -30% more AGN than the Miyaji et al. (2000) sam- 
ple and is comparable to that of Hasinger et al. (2005), including 
not only Type-1 AGN but also sources identified as Type-2 AGN 
in order to gain an insight into the evolutionary properties of the 
whole AGN population at soft X-rays. 

Barger et al. (2005) measured the cosmic evolution of AGN 
in the Chandra Deep Field inthe 2-8 keV energy band finding it 
consistent with a Pure Luminosity Evolution model. Silverman 
et al. ( 120081 studied the XLF of high redshift AGN in the 2- 
8 keV band but they did not consider the intrinsic absorption of 
the sources due to the limited count statistics. Other works in 
the 2-10 keV band such as Ueda et al. d2003) and La Franca 
et al. (120051 ) have taken into account the Nh distribution of the 
sources when calculating the XLF, but an important fraction 
of the Nh values were derived from the hardness ratios of the 
sources rather than from a spectral analysis. This may introduce 
a certain degree of uncertainty since a template spectrum has to 
be assumed to compute the Nh- The XLF of very hard sources 
(>5 keV) is almost unexplored so far with the exception of the 
work by Delia Ceca et al. (fJ008 ), who studied the de-evolved 
(z = 0) luminosity function of absorbed and unabsorbed AGN in 
the 4.5-7.5 keV band but were unable to perform detailed evolu- 
tionary studies of the absorbed population due to the low statis- 
tics. 

In this work we use the XMM-Newton Medium Survey 
(XMS, Barcons et al. 120071 1. along with other highly complete 
deeper and shallower surveys, to compute the X-ray luminosity 
function in several energy bands. Furthermore, given the avail- 
ability of high-quality X-ray spectral information in the XMS, 
we are able to model the intrinsic absorption of the hardest 
sources (4.5-7.5 keV) as a function of the X-ray luminosity up to 
column densities of ~ 10 24 crrT 2 . These issues are key tools to 
probe the accretion history of the Universe across cosmic time. 
Thanks to the extremely high identification completeness of the 
XMS sample (-96% in the 0.5-2 keV band, and -85% in the 
2-10 keV and 4.5-7.5 keV bands) and the accompanying surveys 
we have assembled an overall sample of -1000 identified AGN 
in the 0.5-2 keV, -450 identified AGN in the 2-10 keV band, 
and -120 identified AGN in the 4.5-7.5 keV bands, leading to 
one of the largest and most complete sample up to date in all 
three energy bands. 

This paper is organized as follows: in section [2] we sum- 
marize all the X-ray samples used in this work, giving a short 
overview and details on the level of completeness, sky areas, 
flux limits and energy bands in which they have been selected. 
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Fig. 1. Luminosity-redshift plane of different X-ray surveys in 
the Soft (top panel), Hard (Center panel) and Ultrahard (Bottom 
panel) bands. 



In section[3]we study the absorbing column density distribution 
of the spectroscopically identified AGN detected in hard X-rays 
(>2 keV) and we model the intrinsic fraction of absorbed AGN 
as a function of both luminosity and redshift in the 2-10 and 4.5- 
7.5 keV bands. In section H] we compute the X-ray luminosity 
function of AGN in several energy bands using two methods: a 
modified version of the classic l/V a method to construct binned 
luminosity functions, and a Maximum Likelihood fit to an an- 
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alytic model using all the sources available without binning. In 
section [5] we discuss the results obtained and we compare them 
with those obtained in previous works. Finally, the conclusions 
extracted from this work are reported in section[6] 

In what follows, X-ray luminosities and intrinsic column 
densities Nh are in units of erg s _1 and crrT 2 , respectively. 
Throughout this paper we have assumed a cosmological frame- 
work with H = 70 km s -1 Mpc -1 , Q M = 0.3 and O a = 0.7 
(Spergelatal.|2003). 



2. The X-ray data 

The backbone of the work presented in this paper is the XMS 
sample (Barcons et al. |2007). This survey has typical exposure 
times of about 15 ks and hence is a survey reaching moderate 
depths, but in order to cover a wide luminosity and redshift range 
it is imperative to combine it with other complementary shal- 
lower and deeper surveys. Shallower, wider area surveys will 
provide significant numbers of bright sources at low redshifts, 
while deep pencil-beam surveys will probe fainter sources at 
greater distances. 

Since the XMS survey is composed of different subsamples 
selected in several energy bands, we need to incorporate addi- 
tional surveys for each of them. In the analysis performed in this 
paper we have used three energy bands defined as follows: 

- Soft: 0.5-2 keV 

- Hard: 2-10 keV. 

- Ultrahard: 4.5-7.5 keV 

In this section we will describe all the surveys used in 
this paper, which are also summarized in Table Q] They are 
all well-defined flux-limited surveys with very high identifica- 
tion completenesses, taken from already published catalogues 
of sources from past and present X-ray observatories (XMM- 
Newton, Chandra, ASCA, ROSAT). In order to have an opti- 
mal L x - z plane coverage (see Figure [TJ, we will constrain our 
analysis to those sources identified as AGN with redshifts in the 
range 0.01 < z < 3 in the Soft and Hard bands (spanning up to 
seven and six orders of magnitude in luminosity, respectively), 
and 0.01 < z < 2 in the Ultrahard band (spanning up to four or- 
ders of magnitude in luminosity), being the vast majority of them 
spectroscopical redshifts. The total sky areas covered by these 
surveys are shown in Figure [2] For the purposes of this work, 
and to avoid errors and biases caused by further classification, 
we have used the entire AGN population available (within the 
redshift limits stated above) irrespective of whether they were 
optically identified as Type-1 or Type-2 AGN. 



2. 1 . XMM-Newton Medium Survey 

The XMM-Newton Medium Sensitivity Survey (XMS, Barcons 
et al. I2007I ) is a flux-limited survey of serendipitous X-ray 
sources at intermediate fluxes built from 25 XMM-Newton fields 
of the AXIS sample (Carrera et al. |2007l ). It covers a geometric 
sky area of 3.33 deg 2 and comprises 318 distinct X-ray sources, 
out of which 272 (86%) have been spectroscopically identified 
using a number of ground-based facilities. This fraction, how- 
ever, increases significatively in the Soft band, where the identi- 
fication completeness is 96%. The flux limit in the Soft band is 
1.5x10 erg crrT 2 s _1 , while in the Hard and Ultrahard bands 
are 3.3xl0~ 14 and 6.8xl0~ 15 erg cirT 2 s _1 , respectively. In the 
Soft and Hard bands these limits are well above the sensitivity 
of the data, while in the Ultrahard band it corresponds to the flux 
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bands. 



of the faintest detected source in this sample (see Barcons et al. 
120071 for details). For this work we have excluded the sources 
classified as stars, normal galaxies and clusters of galaxies, end- 
ing up with a sample of 178 AGNs in the Soft band, 120 in the 
Hard band and 57 in the Ultrahard band. A fraction of the ob- 
jects classified as normal galaxies (which are ~4% over the total 
sample in the Soft and Ultrahard bands, and ~5% in the Hard 
band) might be optically elusive AGN given their high intrinsic 
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Table 1. Summary of surveys used in this work, along with their flux limits, sky coverage, total number of sources (identification 
completeness in parenthesis), and the number of identified AGN for each energy band. 



Soft (0.5-2 keV) 



Survey 


Flux limit (erg cm 


" 2 s- 1 ) 


Area (deg 2 ) 


Ntolal 


Nagn 


RBS 


2.5xl0-' 2 




20300 


953 (100%) 


310 


RIXOS8 


8.4xl0~ 14 




4.44 


105 (100%) 


40 


RIXOS3 


3.0xl0~ 14 




15.77 


296 (94%) 


182 


XMS 


1.5xl0" 14 




3.33 


210 (96%) 


178 


UDS 


1.2xl0" 15 




0.36 


94 (95 %f 


73" 


CDF-S 


5.5xl0~ 17 




0.12 


307 (99%)" 


226" 


Hard (2-10 keV) 


Survey 


Flux limit (erg cm" 




Area (deg 2 ) 


Ntolal 


Nagn 


AMSS 


3.0xl0~ 13 




68 


87 (99%) 


79 


XMS 


3.3xl0~ 14 




3.33 


159 (84%) 


120 


CDF-S 


4.5x10-'" 




0.11 


251 (99%)" 


236" 


Ultrahard (4.5-7.5 keV) 


Survey 


Flux limit (erg cm" 




Area (deg 2 ) 


N 10 ,al 


Nagn 


HBSS 


7.0xl0~' 4 




25.17 


67 (97%) 


62 


XMS 


6.8xl0~ 15 




3.33 


70 (86%) 


57 



8 Including photometric redshifts (see text). 



luminosities (logLx ^ 42). In these cases we have adopted a 
conservative position and excluded them from the final sample, 
both in the XMS and the other surveys involved. 

2.2. XMM-Newton Hard Bright Sample 

The XMM-Newton Hard Bright Survey (HBSS) is part of a 
bigger survey project known as XMM-Newton Bright Survey 
(XBS, Delia Ceca et al. 120041 Caccianiga et al. 120081) . In partic- 
ular, the HBSS is a survey of sources at high Galactic latitudes 
detected in the 4.5-7.5 keV (Ultrahard) band down to a flux limit 
of 7xl0~ 14 erg cm" 2 s _1 , with a flat sky coverage of 25.17 deg 2 . 
The sample contains 67 sources. All of these sources but two 
have been spectroscopically identified, as reported in Caccianiga 
et al. (2008), with 62 of them being AGN. Along with the spec- 
troscopic redshifts, information on the intrinsic absorption col- 
umn densities for each source is also provided in Delia Ceca et 
al. (f2008) which will allow us to model the absorption of the 
sources detected in the Ultrahard band (most of them Type-2 
AGN) as explained in Section[3] 

2.3. Chandra Deep Field South 

The Chandr a Deep Field South (CDF-S, Giacconi et al. I200T1 
Rosati et al. [2002) is one of the deepest surveys in the Soft and 
Hard bands carried out so far, with a total exposure time of 1 
Ms. The source samples are widely discussed in Bauer et al. 
(2004). We have used the 346 sources reported in Giacconi et 
al. ([2002 ) for a total of of 307 and 25 1 sources in the Soft and 
Hard bands, respectively, covering ~0. 125 deg 2 and ~0. 108 deg 2 
down to flux limits of 5.5xl0~ 17 and 4.5xl0" 16 erg cm" 2 s _I 
in both energy bands. The optical imaging strategy and optical 
counterparts catalogue can be found in Giacconi et al. ([2002 ), 
while the spectroscopic identifications are taken from Szokoly 
et al. ( 120041 1. In both cases the ground-based facilities used were 
the FORS1 and FORS2 instruments at the VLT. The redshifts of 



the inconclusive identifications and unidentified sources (which 
sum up to ~60% of the total CDF-S sample) have been taken 
from photometric redshifts estimations in Zheng et al (2004). 
These estimations make use of 12-band data in near ultraviolet, 
optical, infrared and X-rays along with sets of power law models 
for Type-1 AGN and various templates of galaxies. They have 
been calculated using two parallel models: HyperZ (Bolzonella, 
Miral les & Pello 120001 ) and the Bayesian model BPZ (Bemtez 
2000). The photometric redshifts were checked against the se- 
cure spectroscopic identifications in Szokoly et al. (120041) and 
the COMBO- 17 survey (Wolf et al. |2Wll20U3ll2004] ) to ensure 
their reliability, matching well that of both surveys. The over- 
all number of sources in Giacconi et al. (2002 ) is 346, of which 
145 have secure spectroscopic redshifts, 4 are unidentified and 
the remaining 197 sources have been photometrically identified. 
Photometric redshifts with flags 0.2, 0.3 and 0.4 in Zheng et al. 
(2004) are considered low-quality identifications since they rely 
only in one of the methods described in the paper (16 sources 
in total, 4.7% of the total CDF-S sample). We have also made 
use of these 16 identifications in this work since their number is 
too low to significantly affect the XLF parameters given the total 
size of the samples used here. Since we need a highly complete 
survey at faint fluxes in order not to degrade the overall sam- 
ple completeness and given the accuracy of the estimations, this 
is the only case in which we are using photometric redshifts in 
this work. The total number of sources classified as AGN in the 
CDF-S that we are using in this paper is 226 in the Soft band and 
236 in the hard band. The Nh column densities have been taken 
from Tozzi et al. (|2006). 

2.4. ASCA Medium Sensitivity Survey 

The A SCA Medium Sensitivity Survey (AMSS, Ueda et al. 
2005 ) is one of the largest high Galactic latitude, broad band, 
X-ray surveys to date. It includes 606 sources in the 2-10 keV 
(Hard) band over a sky area of 278 deg 2 with hard band fluxes 
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Fig. 3. Intrinsic Nh versus X-ray luminosity of XMS sources detected in the Hard {Left panel) and Ultrahard (Right panel) bands. 
Dots represent those sources identified as Type-1 AGN while the Triangles are those classified as Type-2 AGN. The dashed line at 
log Nh = 22 marks the standard separation between absorbed and unabsorbed AGN. 



spanning between 10~ 13 and 10~ n erg cirT 2 s _1 . For this work 
we have used a flux-limited subsample of the AMSS in the north- 
ern sky (AMSSn). It has a total sky coverage of 68 deg 2 , and all 
but one of the 87 hard X-ray sources detected down to a flux 
limit of 3xl(T 13 erg cirT 2 s _1 have been spectroscopically iden- 
tified (79 of them being AGN). The Nh measurements have been 
taken from Table 2 in Akiyama et al. (2003). Further details on 
the imaging and spectroscopic identifications can also be found 
in Akiyama et al. J20031 >. 

2.5. ROSAT International X-ray/Optical Survey 

The ROSAT International X-ray/Optical Survey (RIXOS, 
Mason et al. 120001 1 is a medium sensitivity survey of high 
Galactic latitude X-ray sources. The sample contains 401 
sources divided in two subsamples: 64 ROSAT fields with 296 
sources down to a flux limit of 3xl0~ 14 erg cirT 2 s (known 
as RIXOS3), plus 18 further fields containing 105 sources down 
to a flux limit of 8.4xl0~ 14 erg cirT 2 s" 1 (RIXOS8) in the Soft 
band. RIXOS3 covers 15.77 deg 2 in the sky while RIXOS8 cov- 
ers 4.44 deg 2 , summing a total sky area of ~20.2 deg 2 . The 
RIXOS3 subsample has been spectroscopically identified up to 
a rate of 94% (with total population of 182 AGN), whereas the 
RIXOS8 subsample has been completely identified (adding 40 
AGN to the full sample), which makes them ideal for evolution 
studies in the Soft band. 



2.6. ROSAT Deep Survey - Lockman Hole 

The ROSAT Deep Survey (RDS, Hasinger et al. [T9981 collects 
all the observations performed by ROSAT in the period 1990- 
1997 in the direction of the Lockman Hole, which is one of the 
areas of the sky with a minimum of the Galactic Hydrogen col- 
umn density (Lockman et al. |19861 l. In this paper we use an 
extension of t he RD S, the ROSAT Ultra Deep Survey (UDS, 
Lehmann et al. 2001), which comprises 94 X-ray sources down 
to a flux limit of 1.2xl0~ 15 erg cirT 2 s _1 (~5 times fainter than 
the RDS) in the Soft band. This sample is 90% spectroscopically 
identified, containing 70 AGN (mainly Type-1) plus 3 photomet- 
rically identified Type-2 AGN that have been also included in the 
analysis for a total of 73 AGN. 



2. 7. ROSAT Bright Survey 

The ROSAT Bright Survey (RBS, Fischer et al. [19981 Schwope 
et al. 120001) aimed to identify the brightest ~2000 sources de- 
tected in the ROSAT All Sky Survey (RASS, Voges et al. [19991 
at high Galactic latitudes, excluding the Magellanic Clouds and 
the Virgo cluster, with PSPC count rates above 0.2 s _1 . This 
sample contains bright sources (flux limit ~2.5xl0~ 12 erg cm -2 
s _1 in the Soft band) with redshifts mainly below 1 and a total 
sky coverage of ~20300 deg 2 . The sample has been completely 
identified via spectroscopic observations and includes 3 10 AGN, 
plus a strong population of galaxy clusters which has been ex- 
cluded from this work. 



3. The Nh function of hard sources 

This work aims at calculating the cosmological evolution of the 
X-ray luminosity function of all AGN (both unabsorbed and ab- 
sorbed) within our sample in three energy bands. This should be 
done by calculating the intrinsic (before absorption) luminosity 
and absorption (Nh) function in order to obtain results free from 
selection effects. Since hard X-ray photons are less affected by 
absorption, absorbed AGN are more likely to be detected at en- 
ergies above 2 keV and hence we have applied such method to 
the Hard (2-10 keV) and Ultrahard (4.5-7.5 keV) sources, for 
which we have detailed spectral data (photon index T and intrin- 
sic absorbing column densities Nh, Mateos et al. [2005 ). With 
~25% of the XMS sources in the Hard and Ultrahard bands clas- 
sified as Type-2 AGN, we find that they typically log Nh > 22, 
a value commonly used to separate absorbed AGN from the un- 
absorbed ones (see Figure [3j. An appropiate calculation of the 
X-ray luminosity function of these sources would therefore re- 
quire a previous modelling of the intrinsic absorption in order to 
avoid possible selection effects (see Delia Ceca et al. 2008]for a 
discussion on this subject). On the contrary, the sources detected 
in the Soft band (0.5-2 keV) are so strongly affected by absorp- 
tion that, in this band, we can assume that we are effectively 
sampling mainly unabsorbed AGN. 

The Nh function f(Lx,z;Nn) is a probability distribution 
function for the absorption column density as a function of 
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Fig. 4. Observed N H distributions of the joint XMS/AMSS/CDF-S sample in the Hard band (Left panel) and XMS/HBSS sample 
in the Ultrahard band (Right panel), respectively. The distribution is shown for the full sample and in different luminosity ranges. 
The first bin accounts for both the sources observed in that bin (lower step) and these sources plus the sources with log Nh < 20 or 
uncomputed absorptions (higher step). 



the X-ray luminosity and redshift. It is measured in units of 
(log Nh)~ 1 and is normalized to unity over a defined Nh region: 



^logiV 
JlogNu 



f(L x ,z;NH)d\ogN H = \ 



(1) 



We have chosen the Nh limits to be log Nu ml „ - 20 and 
log NH„ m , = 24, given the intrinsic Nh range spanned by our joint 
samples (see Figure|4|. The observed fraction of absorbed AGN 
(those with log Nh > 22) is measured by the parameter iff, which 
is in general function of both luminosity and redshift. If we com- 
pare the value of if/ in different luminosity ranges we can observe 
that the fraction of absorbed AGN is not constant, decreasing as 
the luminosity increases in both the Hard and Ultrahard bands. 
The fraction of absorbed AGN clearly increases with redshift 
in the case of the Hard band (see Figure 0, mainly due to the 
sources above z > 2 coming from the CDF-S survey. If we pay 
attention only to the points below z ~ 2, the evolution in redshift 
looks milder. On the other hand, there is not significant varia- 
tion in the value of ifr with redshift at a given luminosity in the 
Ultrahard band (see Figure[6]i, which can be partly explained by 
the poor coverage in redshift of this sample. Therefore, we will 
assume that the formal expression of tfr is dependent on both the 
X-ray luminosity and redshift and hence we have formalized its 
parametrization as a linear function of log Lx and z, similarly as 
in La Franca et al. 120051 



i/,(L x ,z) = if/44 [(log I* - 44)/? L + 1] [(z - 0.5)& + 1] 



(2) 



where 1/144 is the fraction of absorbed AGN at log Lx = 44 and 
z = 0.5, and /3l and f3 z are the slopes of the linear dependencies 
on luminosity and redshift, respectively. 

Given equationQ] the normalized Nh function can be written 



f(L x ,z;N H ) 



2 



; 20 < logN H < 22 
;22 < logN H < 24 



(3) 



In order to obtain the best-fit values of the free parameters 
^44, Pl and f3 z we have performed a x 1 fit on the sources of 
the Hard and Ultrahard samples, although in the latter case we 



have fixed /5 Z — to account for the absence of dependence 
on redshift observed in this band. The values thus obtained are 
if>44 = 0A1 + _° Q ™, /3 L = -0.22!^ and /3 Z = 0.57+°;^ for the 
Hard band, and ^44 = 0.22 + 0.04 and yS L = -0.45 + _° f 5 for 
the Ultrahard band (see Table |2). The lcr errors correspond to 
A* 2 = 1. 

We find that our best-fit value at logL2-io = 44 (^44 = 
0.41^^) is in excellent agreement with the ones calculated by 
Ueda et al. ( [20031 (^ 44 = 0. 41+0 03) and La Franca et al. ( [20051 ) 
((^44 — 0.42^^). Hasinger (2008) found a strong linear decrease 
from -80% to -20% in the range l ogLx = 42 - 46. A linear fit 
to the data in Table 5 of Hasinger (2008) yields a best-fit value 
at log L2-10 = 44 of 0.38 + 0.04, and a slope of -0.226 + 0.014, 
also in agreement with our results in the 2-10 keV band within 
the error bars. The decrease in the absorbed AGN fraction at 
high luminosities have been reported by many authors (Ueda et 
al. 120031 La Franca et al. |2005l Akylas et al. I2U061 Delia Ceca et 
al. 120081 . 

La Franca et al. (2005) reported that the absorbed fraction 
is also dependent on the redshift. Redshift dependence is also 
found by Treister & Urry (120061 1 and Ballantyne et al. ( 120061) . 
while Dwelly & Page (2006) did not find dependence on either 
the luminosity or redshift. In a recent paper, Hasinger (2008 ) 
has found a strong decrease in the fraction of absorbed AGN 
with X-ray luminosity and a significant increase of that fraction 
with redshift. The result of Hasinger (2008 ) suggests that the 
evolution of this fraction at a fixed luminosity of logL2-io = 
43.75 is faster than the result obtained by Treister & Urry (2006 ) 
(probably due to the fact that the latter used the Broad Line AGN 
classification only, which tends to overestimate the fraction of 
absorbed AGN at low redshifts) and is consistent with our best 
fit model at a typical redshift of z = 0.5 (\f/(z = 0.5) = 0.37+0.07 
against our prediction of if/(z = 0.5) = 0.39 + 0.08). 

The lack of any dependence on redshift of our Ultrahard 
sample, however, must be handled with caution since it comes 
from two similar XMM-Newton surveys only and spans a lim- 
ited redshift range. As stated in Perola et al. d2004| ), who found 
no luminosity dependence from a sample of 1 17 sources in the 
2-10 keV band, deeper X-ray surveys are needed to take into ac- 
count those sources at higher redshifts and lower luminosities in 
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order to fully investigate the true incidence of absorption. Only 
combining several surveys at different depths, as we have done 
in the Hard band, and therefore covering wider regions in the 
Lx - z plane is possible to unveil the real dependencies. 

In a recent work by Delia Ceca et al. (2008 ) the dependency 
on luminosity of the fraction of absorbed AGN detected in the 
4.5-7.5 keV band is calculated using Nh — 4 x 10 21 cirT 2 as 
the dividing value between obscured and unobscured AGN. This 
corresponds to Ay ~ 2 mag for a given Galactic Ay/Nn ratio of 
5.27 x 10 22 mag cirT 2 which is the best dividing value between 
optical Type-1 AGN and optical T ype-2 AGN in the HBSS sam- 
ple as shown in Caccianiga et al. (2008 ). If we repeat the calcu- 
lations in the Ultrahard band using this Nh value as a dividing 
line we obtain ifr 44 = 0.30+jJ ^ and j3 L = -0.42+°-^ as the best 
fit parameters for the Nh function. They are slightly different 
but consistent within the lcr error bars to those calculated with a 
dividing value of Nh = 10 22 cirT 2 . 

Delia Ceca et al. (2008 ) found a fraction of absorbed AGN 
withL 2 -io > 3 x 10 42 ergs-' of 0.57 + 0.11 which was in excel- 
lent agreement with the results obtained by a variety of SWIFT 
and INTEGRAL surveys (see Table 3 in Delia Ceca et al. l2008l . 
The value predicted by our Ultrahard Nh function, converting 2- 
10 keV luminosities to 4.5-7.5 keV using T = 1.7, is 0.55 + 0.18 
using a dividing value of Nh = 4 x 10 21 cirT 2 , which fully 
agrees with that of Delia Ceca et al. (120081 ). Considering a divid- 
ing value of Nh - 10 22 cm" 2 , our predicted fraction decreases 
to 0.42 ± 0.18 which coincides with the INTEGRAL result of 
Sazonov et al. (2007) (0.42 + 0.09) who used the same dividing 
Nh value between absorbed and unabsorbed AGN. 

4. The X-ray Luminosity Function 

In this section we will calculate the X-ray luminosity function 
(XLF) of our sources in the Soft, Hard and Ultrahard bands. The 
differential XLF measures the number of AGN per unit of co- 
moving volume V and log Lx, and is a function of both luminos- 
ity and redshift: 



d®(L x ,z) _ d z N(L x ,z) 
d log L x dVd log L x 



(4) 



and it is assumed that is a continous function over the luminosity 
and redshift ranges where it is defined. 

As a first approach we will estimate the XLF in fixed lumi- 
nosity and redshift bins (section |4"TT > which will allow us to have 
a general overview of the overall XLF behavior. In section 14.21 
we will express the XLF in terms of an analytical formula over 
which we will perform a Maximum Likelihood (ML) fit, using 
the full information available from each single source and thus 
avoiding biases coming from finite bin widths and from the im- 
perfect sampling of the Lx - z plane. 

4.1. Binned luminosity function 

There are a variety of methods to estimate the binned XLF. The 
classical approach is the l/V„ method (Schmidt I19681 ). which 
was later generalized by Avni & Bahcall ( 1980 ) for samples with 
multiple flux limits. This method has been widely used to com- 
pute the binned XLF in evolution studies of flux-limited sam- 
ple s (i.e. Maccacaro et al. [1991] Delia Ceca et al. [T9921 Ellis et 
al. 119961 ) but it can lead to systematic errors, especially when 
dealing with sources very close to the flux limit. A number of 
techniques have been proposed to solve this problem (Page & 
Carrera 2000 Miyaji et al. 1200 11 ). In this work, we have used the 



alternative method proposed by Page & Carrera (2000) to cal- 
culate the differential binned luminosity function, which avoids 
most of the biases that accompany the classic 1 /V a method. 
The differential XLF can be approximated by: 



d®(L x ,z) 
d\ogL x 



N 



m f z dzd\ogL x 



(5) 



where N is the number of objects found in a given luminosity- 
redshift bin ALAz which is surveyed by the double integral in the 
denominator. The differential comoving volume has been calcu- 
lated using the expression from Hogg (119991 ): 



dV _ c_ 
~dz ~H< 



' / D L \ 2 1 , \-l/2 

r 7^- (fl* 1 + z) 3 + n A Q(L X , z) (6) 



where L\ is the luminosity distance and Q(Lx,z) is the solid 
angle subtended by the survey. Note that some portion of the 
ALAz bin may represent objects that are fainter than the survey 
flux limit. The upper limit of the integral in redshift z max (L) is 
hence either the top of the redshift shell Az or the redshift at the 
survey flux limit for a given luminosity: z max (L) = Zu m (L x , S /,„,). 

If the sample is composed of multiple flux-limited surveys 
(as it happens in this work), they are added 'coherently' as ex- 
plained in Avni & Bahcall (1980). This means we will assume 
that each object could be found in any of the survey areas for 
which is brighter than the corresponding flux limit. 

The binned XLF thus obtained for our sample presents an 
overall double power-law shape, with a steeper slope at brighter 
luminosities beyond a given break luminosity that takes place 
somewhere in the range log Lx — 43 - 44. The position of this 
break seems to change with redshift, moving to higher values of 
log Lx as we go deeper, thus showing clear evidences of evolu- 
tion in the AGN that compose our sample. The shape at higher 
redshifts is less constrained due to the limited statistics. An an- 
alytical model is therefore needed to account for these changes 
in shape as it evolves with z, and whose parameters could be 
compared with previous works. 

4.2. Analytical model 

In the light of the results obtained in section l4"Tl we will express 
the XLF as an analytic function with an smooth behavior over 
the range of Lx and z under study in this work. The parameters 
of this function will describe its overall shape and how it evolves 
with luminosity and redshift, which will express physical prop- 
erties of the population under study. 

As shown in many works in the Soft and Hard bands (Miyaji 
et al. 120001 Ueda et al. 120031 Hasinger et al. 120051 Silverman et 
al. 120081 ), the XLF seems to be best described as a double power 
law modified by a factor for evolution. Here we implement two 
different models: the simpler Pure Luminosity Evolution (PLE) 
model in which the X-ray luminosity evolves with redshift, and 
the more complicated Luminosity-dependent Density Evolution 
(LDDE) model in which the evolution factor not only depends 
on the redshift but also on the luminosity. Both models can be 
expressed respectively as: 



d®(L x ,z) d®(L x /e(z),0) 



d log L x 



dlogLx 



for the PLE model, and: 

d®(Lx,z) _ d<t>(L x ,0) 
dlogLx dlogLx 



e(z,L x ) 



(7) 



(8) 
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Fig. 5. Left panel: Fraction of absorbed AGN as a function of the 2-10 keV luminosity. Right panel: Fraction of absorbed AGN as a 
function of redshift in the luminosity range log Lx = 43.5 — 44.5. Dashed lines represent the best-fit model of the Nh function using 



Nh = 10 22 cm 2 as the dividing value between obscured and unobscured AGN. 
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Fig. 6. Left panel: Fraction of absorbed AGN as a function of the 4.5-7.5 keV luminosity. Right panel: Fraction of absorbed AGN 
as a function of redshift in the luminosity range logL^ = 43 — 43.75. Dashed lines represent the best-fit model of the Nh function 



using Nh = 10 22 cm 2 as the dividing value between obscured and unobscured AGN. 



for the LDDE model. The first factor of the second term in equa- 
tions [7] and [8] describes the shape of the present-day XLF for 
which we adopt a smoothly connected double power law: 



d®(L x ,0) 
dlogLx 



= A 



Loj \Lo) 



(9) 



where j\ and are the slopes, Lq is the value of the luminosity 
where the change of slope occurs, and A is the normalization. 
The evolution factor of the PLE model is expressed as: 



;z < z c 



f(l+z) Pl 
e(z) = W)(r^r 



(10) 



where the parameters p\ and p2 account for the evolution be- 
low and above, respectively, the cut-off redshift z c . In the LDDE 
model we will assume that the cut-off redshift depends on the 
X-ray luminosity: 



, r , [(l+z)" 1 ;z<Zc(L x ) 
e(z, L x ) = < , . ( i +z \P2 , 



(11) 



Zc(Lx) 



;L x >L a 
;L x <L a 



(12) 



where a measures the strength of the dependence of z c with 
luminosity. 

4.2.1 . Model fitting to Soft sources 

We have fitted our sample in the Soft band to the PLE and 
LDDE models described above using a Maximum Likelihood 
(ML) method, which optimally exploits the information from 
each source without binning and is therefore free from all the 
biases commented in section |4~T1 In this band we are sampling 
mainly unabsorbed AGN and therefore the derived XLF is that 
of the unabsorbed population in the 0.5-2 keV band. 

The ML technique we have implemented here is that of 
Marshall et al. ( I1983I ), in which the likelihood function is de- 
fined as the product of the probabilities of observing exactly one 
object in the differential element dzd log Lx at each (zi,L Xi ) for 
the jV objects in the sample, and of the probabilities of observ- 
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Soft (0.5-2 keV) Hard (2-10 keV) 


Ultrahard (4.5-7.5 keV) 




PLE 


LDDE PLE LDDE 


PLE 


LDDE 


Present day XLF parameters 


A" 

logL 6 
7i 


50.73±3.79 

42 90 +0 03 

_ 002 

0.34±0.01 


3.76±0.38 17.96 + ^ 4.78^0 
43.56±0.05 43.60+0.13 43.91^ 
0.72±0.02 0.81 ±0.06 0.96+0.02 


10.41 + !^ 
43.00 + |} 2 f 
1.14+0.01 


1.32+0.20 

43 43+°- 23 

i 90+O.O8 


72 


2.01 ±0.04 


2.04±0.04 2.35±0.07 


9 79+0.29 
L - -0.28 


2.53+0.32 


Evolution parameters 


Pi 

Pi 


1.78±0.06 
0.00 (fixed) 


3.38±0.09 2.04 + °| 3 4.07 + °$ 
-1.5 (/aerf) 0.00 (fixed) -1.5 (/terf) 


3 03 +0 ' 17 

-0 22 

0.00 (/ted) 


6 46 +0 69 

-0.29 

-1.5 (fixed) 


Zc 


1.7 (fixed) 


1 .42 (/weed) 1 .9 (//xerf) 1 .9 (/ixerf) 


1.9 (fixed) 


1.9 (/ixcrf) 


logL fl * 




44.6 (/wed) ... 44.6 (fixed) 




44.6 (fixed) 


a 




0.100+0.005 ... 0.245±0.003 




0.245 (/wed) 


N H function parameters 


Pl 

A 




041 +0.03 

o.ti _ 004 
n 99+O.04 

57 +( ™ 

—0.10 


0.22+0.04 
-0 45+ - 20 

"'^ -0.25 

0.00 (fixed) 


P2DKS(L X ,ZY 


2xl0~ 4 


0.21 3xl0~ 4 0.18 


0.28 


0.84 



a In units of lO -6 ^ Mpc" 3 . 
b In units of fc? erg s -1 . 
c 2D K-S test probability. 



ing zero objects in all other differential elements in the accessible 
regions of the z- Lx plane. The expression to be minimized is 
hence: 



J — — ^ Zj,-i "< — n — ; — 

-2 zSr r £ 2 ^^(L X , z) ^)^iogL x 

Here, the index / runs over all the sources present in the sam- 
ple whereas the index j runs over all the surveys that compose 
the sample. In case there is incompleteness in the spectroscopic 
identifications in the j-th survey, the factor Ci(Lx,z) accounts 
for the completeness of the identifications at a given X-ray flux. 
Here we assume that the redshift distribution of the unidentified 
sources is the same as that of the identified sources at similar 
fluxes and hence the effective survey area is the geometric area 
multiplied by this factor. Note that this assumption is not correct 
when the source is unidentified by non-random effects, but we do 
not expect this to affect our results significatively given the high 
degree of completeness achieved by our combined sample. The 
integrals are calculated over the full redshift (0.01 < z < 3) and 
luminosity (40 < log Lx < 46) ranges spanned by our sample in 
the Soft band taking into account the flux limits of the different 
surveys that compose the sample. There are 5 sources with lumi- 
nosities in the range (40 < log Lx < 42) that have been classified 
as AGN (see Szokoly et al. 2004). Although they could be nor- 
mal galaxies misclassified as AGN, we have decided to trust the 
original classification and include them in the final sample. The 
remaining sources labeled as AGN in this range come from the 
photometric classification in Zheng et al. (2004). 

The expression for S is minimized using the MINUIT soft- 
ware package (James |1994| l from the CERN Program Library, 
lcr errors for each parameter are calculated by fixing the param- 
eter of interest at different values and leaving the other parame- 
ters to float freely until AS — 1 . Since the model is rather complex 



and some of the parameters are unconstrained when performing 
the fit using the 9 free parameters, we have fixed some of them. 
In particular we have fixed p2, z c and L a to the values obtained 
by Hasinger et al. (2005) in the 0.5-2 keV band. This left us with 
6 free parameters to fix: A, j\, j2, Lq, p\ and a (see Tabled). 

4.2.2. Model fitting to Hard and Ultrahard sources 

For the sources in the Hard and Ultrahard bands we will make 
use of the information obtained in section [3] when performing 
the ML fit. Given the coverage of the three-dimensional space 
Lx - z - Nh spanned by our Hard and Ultrahard samples, we 
were not able to fit simultaneously the Nh function and the XLF 
(which would lead us to a functional form with 1 1 free param- 
eters). Instead, we will add the Nh function to the analytic ex- 
pression to be minimized but fixing their parameters, i/r 44 , Pl and 
[i z , to those obtained in section [3] Taking this into account, the 
expression for S is: 

S--^l l ln d -^f(L i x ,z'-,N H ) 

^££Qf(L x ,z;N H )^ (14) 
xO(L x , z, N H ) dVHL ^' NH) dzd log L x d log N H 

Again, the integrals are calculated over the full Lx - z - Nh 
space in the ranges 0.01 < z < 2, 20 < logN H < 24, 40 < 
logLx(2-iO) 5 46 (in the Hard band) and 42 < log Lx(4.5-7.5) < 
45 (in the Ultrahard band). The factor C^(Lx,z,Nh) accounts 
for the incompleteness of the identifications and Nh measure- 
ments. Similarly as in the Soft band, there are 2 sources with 
log 

Lx(2-W) < 41 that were originally classified as AGN which 
we have decided to keep in the final sample. As we did in the 
previous section, we fixed the XLF parameters p2, z c and L a to 
those obtained by Ueda et al. ( 120031) for sources detected in the 
Hard band and, in addition, given the poor coverage of the L x - z 
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Fig. 7. X-ray luminosity function in the Soft band in the redshift ranges 0.01-0.5 (Top left panel), 0.5-1.0 (Top right panel), 1.0-2.0 
(Bottom left panel) and 2.0-3.0 (Bottom right panel). The points are the results from the binned XLF and the lines are the best- 
fit curves to our analytical LDDE (solid line) and PLE (dashed line) models evaluated at the median of each redshift range. The 
dot-dashed line is our LDDE model evaluated at z — 0. 



plane achieved by the Ultrahard sample, we have also fixed the 
strength of the dependence of z c on luminosity a to the value we 
obtained when fitting the Hard sample. 

5. Discussion of the results 

5.1. Comparison with other works 

The best-fit parameters of the XLF calculated by the ML method 
explained above are summarized in Table|2] Best fit results to the 
PLE and LDDE models along with data points from the binned 
XLF are shown in Figure s[7|H and [9] 

At first sight our best-fit models in the Soft band (0.5-2 keV) 
reveal strong differences between the PLE and LDDE models. 
Although the PLE best-fit seems to be a fair description of the 
XLF at low redshifts (z < 1), it clearly fails to reproduce the be- 
haviour at high redshifts, where the LDDE model matches better 
with our binned XLF data points. 

The best-fit parameters obtained by our ML technique 
slightly differ from the soft XLF models reported in Miyaji et 
al. d2D00b and Hasinger et al. d2005l > using a variety of ROSAT, 
XMM-Newton and Chandra surveys. For instance, our best-fit 
PLE model requires less evolution (p\ parameter) than that of 
Hasinger et al. (120051 ) model and the present-day luminosity 



break Lq is displaced to fainter luminosities in our model (but 
consistent within 2<x with that of Hasinger et al.). In the Hard 
band, our best-fit PLE model is consistent within the lcr er- 
ror bars with the results obtained by Barger et al. (|2005) in the 
Chandra Deep Field. 

It has been shown in different works that the LDDE model 
for the XLF provides the best framework that describes the evo- 
lutionary properties of AGN, both in soft (Miyaji et al. 2000 
Hasinger et al. 2005) and hard X-rays (Ueda et al. 2003] La 
Franca et al. 120051 1, as well as in the optical range (Bongiorno 
et al. l2007l . 

The LDDE model provides a better overall description of the 
XLF. However, our present day XLF power law shape has flatter 
slopes and the luminosity break Lo is also displaced to fainter lu- 
minosities (similarly as in the PLE model) than in Hasinger et al. 
2005 by more than 2<x, although it is consistent with the LDDE 
models of Miyaji et al. (2000) within the error bars. Our best- 
fit evolution parameter p\ reveals less evolution below the cut- 
off redshift (fixed to the value obtained by Hasinger et al. 2005 
z c — 1 -42) and a much weaker dependence of z c on luminosity a 
than in these works. The significative less evolution p\ found in 
this work with respect that of Hasinger et al. ( 2005 ) is somewhat 
surprising since both samples have comparable sizes. The differ- 
ences are significative even if we repeat the fits putting the fixed 
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Fig. 8. X-ray luminosity function in the Hard band in the redshift ranges 0.01-0.5 (Top left panel), 0.5-1.0 {Top right panel), 1.0-2.0 
(Bottom left panel) and 2.0-3.0 {Bottom right panel). The points are the results from the binned XLF and the lines are the best- 
fit curves to our analytical LDDE (solid line) and PLE (dashed line) models evaluated at the median of each redshift range. The 
dot-dashed line is our LDDE model evaluated at z = 0. 



parameters from Hasinger et al. (2005 ) at their +lcr values, thus 
obtaining p\ best fits ranging from ~3.4 to ~3.9. However, we 
must note that the slopes p\ and p2 in the Hasinger et al. LDDE 
model are dependent on the X-ray luminosity while ours are not. 
For instance in the log Lx = 42 - 46 interval, p ; in the Hasinger 
et al. model ranges between ~3.3 and ~6. 1 . Since our XLF is ex- 
tended towards lower luminosities than that of Hasinger et al., an 
overall lower p\ parameter is consistent with a X-ray luminosity 
dependency. 

Similarly, in the Hard band (2-10 keV) the LDDE model out- 
performs the PLE model at z > 0.5, failing the latter to describe 
the evolution of the binned data points at the faint end of the 
XLF underestimating the data at log Lx > 44 by a factor of sev- 
eral. Comparing our results in the Hard band (2-10 keV) with 
previous major works in this band (Ueda et al. 120031 La Franca 
et al. 120051 Silverman et al. 2008 ) we find an overall good agree- 
ment between our work and theirs, although with few differ- 
ences. The general shape of the present day XLF, described by 
the smooth double power law, is similar to that reported in other 
works within the error bars. The evolution below the redshift 
cut-off p\ is also consistent with the results of Ueda et al. (2003 ) 
and La Franca et al. (2005 ) (both using samples corrected by the 
intrinsic absorption of their sources), and Silverman et al. (2008 ) 
well within the l<x confidence level. The stronger difference be- 



tween these models arises with the strength of the dependence of 
the cut-off redshift z c on luminosity, measured by the parameter 
a. The values of a obtained by Ueda et al. (2003 ) and Silverman 
et al. (12.0081 ) are consistent with each other while the ones cal- 
culated by La Franca et al ((2005 ) (a ~ 0.2) and in this work 
(a ~ 0.25) are sistematically lower. Overall, our best-fit LDDE 
parameters are more constrained than in the other works accord- 
ing to the computed \cr error bars. We have left the evolution 
parameter above the cut-off redshift p2 fixed to that of Ueda et 
al. (the value obtained in La Franca et al. reveals less evolution 
beyond z c but consistent with Ueda et al. within l<x). The model 
in Silverman et al. (2008 ) requires a much stronger evolution 
P2 = -3.27+Q34 than the others. It must be noted, nevertheless, 
that it is extremely difficult to properly constrain the faint end 
of the luminosity function given the necessity of highly com- 
plete deep pencil-beam surveys that account for the population 
of high-redshift low-luminosity AGN. 

In a recent work, Delia Ceca et al. (2008) have computed 
the intrinsic present day XLF of the HBSS sample in the 4.5- 
7.5 keV (Ultrahard) band for absorbed and unabsorbed AGN. 
In this paper we have added the HBSS sample to ours, and per- 
formed the same analysis as in the Soft and Hard bands using 
both absorbed and unabsorbed sources. This way we have dou- 
bled the number of available AGN and improved the Lx~z cov- 




Fig. 9. X-ray luminosity function in the Ultrahard band in the redshift ranges 0.01-0.5 (Top left panel), 0.5-1.0 (Top right panel) and 
1.0-2.0 (Bottom panel). The points are the results from the binned XLF and the lines are the best-fit curves to our analytical LDDE 
(solid line) and PLE (dashed line) models evaluated at the median of each redshift range. The dot-dashed line is our LDDE model 
evaluated at z — 0. 



erage, which is directly reflected in that the best-fit parameters 
are more constrained and the associated error bars are signifi- 
cantly smaller. Moreover, given the availability of detailed spec- 
tral information for the vast majority of the sources in both sam- 
ples, we have carried out the XLF analysis by convolving the 
Nh function with the XLF analytical model when performing 
the fit (see section |4.2.2t . Like in the other energy bands, the 
PLE best-fit clearly underestimates the data at faint luminosities 
and high redshifts. From the binned data points it can be inferred 
that a very strong evolution p\ is required to describe properly 
the behaviour of the Ultrahard sources, which is achieved by 
the LDDE model. Our results are fully consistent with those re- 
ported in Delia Ceca et al. (2008) albeit with smaller error bars. 
The shape of our best-fit present day XLF would correspond to 
the average of the absorbed and unabsorbed present day XLF 
calculated in Delia Ceca et al. thus lying inside the lcr confiden- 
dence levels spanned by their parameters. Our best-fit value for 
the evolution parameter p\ — 6.46^29 is in excellent agreement 
with that of Delia Ceca et al. (2008 ) (p\ = 6.5) and also with that 
of Bongiorno et al. (2007 ) obtained from a selected AGN sample 
in the optical range from the VIMOS-VLT survey (p\ = 6.54) 
and that of the bolometric quasar luminosity function of Hopkins 
et al. (120071 1 (p\ — 5.95 + 0.23). This value represents a much 
stronger evolution below the cut-off redshift than that obtained 



by Ueda et al. (2003) (p x = 4.23 + 0.39) and La Franca et al. 
(120051 1 (pi = 4.62 + 0.26) in the 2-10 keV band. 

5.2. Accretion history of the Universe 

The X-ray luminosity function of AGN, whose bolometric emis- 
sion is directly linked to accretion power, constrains the history 
of the formation of the supermassive black holes that reside in 
the galactic centers along cosmic times. 

The comoving density of AGN as a function of redshift 
can be calculated straightforward from our best-fit XLF LDDE 
model: 

(15) 

J L] dlogLx 

As we can see in Figure [TO] the most luminous AGN are 
formed before the less luminous ones in all bands. AGN with 
log Lx > 44 reach a maximum in density at redshift ~ 1 .5 while 
fainter AGN (\ogL x < 44) peak at z ~ 0.7. 

Similarly to the comoving density, we can derive the lumi- 
nosity density as a function of redshift in all bands and, from 
that, calculate the accretion rate density. Here we will assume 
that the accretion rate onto a supermassive black hole is related 
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Fig. 10. Comoving density of AGN in different luminosity bins 
as a function of redshift in the Soft (top panel), Hard (Center 
panel) and Ultrahard (Bottom panel) bands. Overplotted are the 
predictions from our best-fit LDDE model. 



to the bolometric luminosity by a constant factor e which is the 
radiative efficiency of the accretion flow (Marconi et al. 2004 
La Franca et al. 120051 : 



Lboi = 



1 -e 



-M acc c z 



(16) 



We can derive the bolometric luminosities by means of a 
bolometric correction factor K simply using L/, G / — KL X - The ac- 
cretion rate density is hence (Soltan [T9821 Marconi et al. 2004): 



Paced) 



i-£ r Li 



KL. 



d<&(L x ,z) 
d\ogL x 



dlogL x 



(17) 



We will assume a nominal radiative efficiency of e = 0.1 (Yu 
& Tremaine 120021 Marconi et al. |2004l Barger et al. 120051 ) . The 
values of the bolometric correction K are derived from the poly- 
nomic expressions of Marconi et al. (2004), that accounts for 
changes in the overall spectral energy distribution of AGN as a 
function of the optical luminosity. Hence, the total accreted mass 
onto supermassive black holes is: 



P(z) 



-r 



c (z)^-dz 
dz 



(18) 



where we assume that the initial mass of seed black holes at zo 
is negligible with respect to the total accreted mass (La Franca 
et al. EoTET l. 

In Figures QT| and [12] is clearly seen that the vast majority of 
the accretion rate density and the total mass density are produced 
by low-luminosity AGN (logLx < 44). The total accreted mass 
density at z = obtained from our XLF model (~ 3 x 10 5 M Q 
Mpc~ 3 ) in the hard band is in good agreement with that of La 
Franca et al. (l2005l) (3. 2hl {) x 10 5 M Q Mpc" 3 ) derived from their 
intrinsic XLF, or the local black hole mass function of Marconi 
et al. d20Q4| > (4.6±j^ x 10 5 M Mpc~ 3 ) determined applying 
the correlations between black hole mass, bulge luminosity and 
stellar luminosity dispersion to galaxy luminosity and velocity 
functions in the Local Universe. 

From these results it can be inferred that high-luminosity 
AGN grow and feed very efficiently in the early Universe and 
are fully formed at redshifts 1.5-2, whereas the low-luminosity 
AGN keep forming down to z ~ 1. These results are hence in 
full agreement with an anti-hierarchical black hole growth sce- 
nario as shown by the LDDE model of the XLF in a number of 
previous works (e.g. Ueda et al. 120031 Merloni 2004 Marconi et 
al. 120041 La Franca et al. |2005l Hasinger et al. |2005l ). 

6. Conclusions 

We have discussed here the cosmic evolution of a sample of 
AGN in three X-ray bands: Soft (0.5-2 keV), Hard (2-10 keV) 
and Ultrahard (4.5-7.5 keV). The backbone of our sample is the 
XMS survey (Barcons et al. |2007l > which is a flux-limited highly- 
complete sample at medium fluxes. We have combined the XMS 
with other shallower and deeper highly complete X-ray surveys 
in all three bands to end up with a total sample of ~ 1000, 435 and 
119 AGN in the Soft, Hard and Ultrahard bands, respectively. 

We have used the spectral information on the sources that 
compose the Hard (2-10 keV) band to model their instrinsic ab- 
sorption (Nh function). We find dependency of the fraction of 
absorbed AGN on both the X-ray luminosity and redshift. Our 
predictions on the behaviour of the fraction of absorbed AGN in 
this band is in excellent agreement with the results of Ueda et 
al. ( 120031 ) and La Franca et al. ( 12005I ). In recent work, Hasinger 
( 2008 ) studies the evolution of the absorption properties of AGN 
using a compiled sample of 1290 AGN in the 2-10 keV band, be- 
ing our results in good agreement with his within lcr. The same 
analysis was applied to the Ultrahard samples, XMS and HBSS 
(Delia Ceca et al. [2008 ). We find that the fraction of absorbed 
of AGN in the 4.5-7.5 keV band, assuming a dividing value of 
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Fig. 11. Accretion rate density of matter onto supermassive black Fig. 12. Total accreted mass onto supermasive black holes as a 
holes as a function of redshift in the Soft (top panel), Hard function of redshift in the Soft (top panel), Hard [Center panel) 



(Center panel) and Ultrahard (Bottom panel) bands. 



and Ultrahard (Bottom panel) bands. 



Nh = 10 22 cm~ 2 , is dependent on the X-ray luminosity but not 
on the redshift. This could be motivated by the narrow redshift 
range spanned by the sample, with the bulk of the sources lo- 
cated at low redshifts (z < 1). In the same band, Delia Ceca et al. 
(2008 ) suggest a possible evolution of the absorbed AGN popu- 
lation with redshift by comparing their results with other survey 
projects. 



We have calculated the X-ray luminosity function of AGN 
using two methods. First, a modified version of the 1 /V a method 
(Schmidt 1968 ) discussed in Page & Carrera (2000) to com- 
pute the binned XLF. Secondly, a fit to an analytic model us- 
ing a Maximum Likelihood technique (Marshall et al. 1983) that 
fully exploits the available information on each individual source 
without binning. The adopted model consists of a smoothly 
connected double power law that accounts for the present day 
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(z — 0) XLF modified by an evolution factor that depends on the 
redshift (PLE model) or on both redshift and luminosity (LDDE 
model). We have found that the LDDE model outperforms the 
PLE model at faint luminosities and high redshifts. 

In general, our computed XLF in the Soft, Hard and 
Ultrahard bands are in good agreement with the results of other 
published surveys. The best-fit XLF parameters in the Soft sam- 
ple show slight discrepancies in both the overall shape and evo- 
lution with respect previous works in that band (Miyaji et al. 
2000 Hasinger et al. 120051 1. In the Hard band, where we have 
computed the intrinsic XLF taking into account the intrinsic ab- 
sorption Nh of each source, we are in good agreement with the 
results of La Franca et al. ( 120051 1 although there are differences 
in detail with Ueda et al. (120031 1 or Silverman et al. ((2008). Our 
best-fit model shows weaker evolution of the AGN below the 
cut-off redshift than in these works albeit with smaller error bars. 
In the Ultrahard band, we have also calculated the intrinsic XLF 
finding similar results as in Delia Ceca et al. (2008) but with our 
best-fit parameters much more constrained. The results in this 
band show that Ultrahard AGN present a significantly stronger 
evolution below the cut-off redshift than those detected at softer 
energies. A possible explanation for this could reside in the ex- 
istence of a dependency on the X-ray luminosity for the evolu- 
tionary parameter p\ that we have not taken into account in our 
model due to the low statistics in this band. Such a dependency 
has been found in the XLF analysis in other bands (e.g. Ueda et 
al. 120031 Hasinger et al. 2005) showing that p\ linearly increases 
with luminosity. Given that in the Ultrahard band we lack of deep 
surveys, which account mainly for intrinsically faint objects, our 
result could be biased towards high-luminosity sources hence 
providing higher evolution rates below the cut-off redshift. In all 
three bands, the high-luminosity AGN (log Lx > 44) are formed 
before than the low-luminosity ones (log Lx < 44), reaching the 
former a maximum in density at redshift z ~ 1 .5 whereas the 
comoving density of the latter peak at z ~ 0.7. 

Finally, we have used our best-fit XLF to compute the accre- 
tion rate density and total accreted mass onto supermassive black 
holes as a function of redshift. The total black hole mass density 
at z — predicted by our best-fit model is in agreement with that 
computed by La Franca et al. (2005) and the local supermassive 
black hole density derived by Marconi et al. (12004). Although 
the value derived by Marconi et al. (|2004| l is slightly higher than 
those computed by La Franca et al. ((2005 ) and in this work, it 
has to be taken into account that hard surveys like the ones used 
here still miss the Compton-thick AGN population which are 
contributors to the relic black hole mass function. Marconi et al. 
(2004) estimate the contribution of Compton-thick AGN to the 
observed BH mass function in a factor of ~ 1 .5 which would fully 
explain this discrepancy. Similarly, Delia Ceca et al. (2008) esti- 
mate a density ratio between Compton-thick and Compton-thin 
AGN of 1.08+0.44 at logL x ~ 43. 

As predicted by the XLF LDDE model, the high-luminosity 
AGN have a more efficient growth in the early stages of the 
Universe and are fully formed at z ~ 1 .5 - 2 while the less lu- 
minous AGN keep forming down to redshifts below 1 (see e.g. 
Merloni [20041 Marconi et al. 120041 Hasinger et al. I2005l among 
others). This behaviour is found in all energy bands under study, 
thus confirming that the evolution of the intrinsic XLF along cos- 
mic time is not caused by changes in the absorption environment 
but by intrinsic variations in the accretion rate at different epochs 
of the Universe. 
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